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ABSTRACT 


1.  INTRODUCTION 

/ 

This  paper  presents  a  numerical  method  for  obtaining  the  super¬ 
sonic,  laminar  vise'  .s  flow  about  arbitrary  geometries  without 
compression  surfaces  at  high  angle  of  attack.  In  particular, 
results  are  presented  for  blunt  biconic  bodies  with  windward 
and  leeward  cuts.  The  basic  approach  used  is  to  solve  the 
steady  three-dimensional  "Parabolized  Navier-Stokes  Equations" 
(PNS) ,  first  derived  for  circular  cones  by  Lubard  and  Helli- 
well^~>2Ltsee  Figu xe-J&f.  These  equations  have  been  used  to 
predict  the  flowfield  for  a  variety  of  different  problems,  in¬ 
cluding  flow  over  sharp  and  blunt  cones  at  angle  of  attack  up 
to  40°  flow  over  spinning  cones  at  angle  of  attack  ^r, 

and  flow  over  cones  with  mass  transfer  and  temperature  varia¬ 
tion  at  the  surface.  In  addition  to  these  results  which  were 
confined  to  circular  cones,  some  limited  results  have  been 
obtained  for  biconic  geometries  ,  non-circular  cones'^  and 
the  NASA  space  shuttle  ^  . 


For  the  three  dimensional  geometries  of  interest  (e.g.  Figure 
1) ,  we  have  developed  a  non-orthogonal  curvilinear  coordinate 
system  that  appears  ideal  for  this  problem.  The  governing 
equations  are  written  in  terms  of  the  non-orthogonal  coordi¬ 
nates  and  their  metrics.  The  three  velocity  components  are 
also  defined  in  the  non-orthogonal  coordinate  directions. 

This  is  different  from  writing  the  equations  in  an  orthogonal 
coordinate  system  and  explicitly  performing  a  coordinate 
trans  forma  t ion . 

Predictions  using  the  new  approach  have  been  obtained  for 
the  body  illustrated  in  Figure  1  at  angles  of  attack  up  to  10° 
and  a  Mach  number  of  10.  Surface  pressure  and  heat  transfer 
have  been  compared  with  wind  tunnel  data  and  show  very  good 
agreement . 
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2.  COORDINATE  SYSTEM  AND  EQUATIONS 


The  steady  Navier-Stokes  Equations  are  first  written  in 
general  curvilinear  coordinates ,  and  then  a  par  icular  co¬ 
ordinate  system  is  chosen.  One  of  the  coordinates,  is 
chosen  in  the  general  axial  direction,  another  C2»  in  a 
direction  away  from  the  body,  and  the  third,  $3 ,  around  the 
body  (Figure  2)  .  It  is  assumed  that  coordinates  with  these 

characteristics  and  that  the  types  of  flows  to  be  studied 

d  3  9 

permit  the  approximations  g^Y<<:  viscous  terms, 

The  resulting  system  is  called  the  "parolized  Navier-Stokes" 
(PNS)  equations . 


In  order  to  achieve  validity  of  the  PNS  approximation,  good 
finite  difference  approximations  to  the  derivatives,  and 
numerically  accurate  boundary  conditions,  it  is  desirable  to 
have  a  coordinate  system  that  has  both  body  and  shock  surfaces 
as  coordinate  surfaces  along  with  the  £2  coordinate  orthogonal 
to  these  surfaces.  To  generate  an  orthogonal  three-dimensional 
coordinate  system  with  these  properties  would  be  a  formidable 
task.  If  it  is  permitted  that  the  5^  and  £3  coordinates  not 
necessarily  be  orthogonal,  but  are  orthogonal  to  the  £2  co¬ 
ordinate,  then  a  coordinate  system  with  the  above  properties 
can  be  generated.  In  summary  we  construct  the  £2*  and  £3 
coordinates  so  as  to  have  the  following  properties: 


1) 

2) 

3) 

4) 


The  body  is  a  coordinate  surface  (52  =  0) . 

The  and  £3  coordinates  are  necessarily  orthogonal 
only  at  the  body. 

The  C2  coordinate  is  orthogonal  to  and  £3* 


The  bow  shock  is  a 


£3  coordinate  surface 
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In  order  to  generate  these  non-orthogonal  coordinates,  we 
start  at  the  body  and  calculate  out  to  the  shock.  The  surface 
of  the  body  is  written  in  cylindrical  coordinates  as 

r  =  r  (z,4>)  (1) 

and  is  defined  to  be  a  vector  in  the  direction  with  the 

coordinate  metrics  g.  .  =  g.  .  g,.  The  £.  and  £_  coordinates 

X  ]  X  J 

on  the  body  (£2  "  0)  are  generated  by  using  the  method  out¬ 
lined  by  Blottner  and  Ellis  in  Reference  (10)  . 

The  coordinates  in  the  region  from  the  body  surface  to  the 
shock  surface  are  obtained  by  taking  the  cross  product, 

23  x  3i*  '  A'*  — '  are  un^t  vectors  in  cylindrical  co¬ 

ordinates,  then  the  ^coordinate  vectors  are  given  by 

3m  =  %T~  *•'  +  ■§?”  j*  +  r  k'  ,  m  «  1,2,3.  (2) 

vm  m  ^m 

Defining  n  as  the  cross  product  between  and  c^,  and  using 
Equation  (2) ,  we  obtain 


where 


Now  set 


n  =  n^  i'  +  nj  i'  +  n^  k '  I 

_  „/3 r  3 4>  34i  3r  \ 

"1  -  ^  51^  -  5I7  m)  • 

_  ( 34>  3z  3z  3<t>  \ 

"2  -  nq  -  jq  Ttq)  -  ' 

„  /3z  3r  3r  3z  \ 

"3  ’  7^  TC[)  • 

^  and  again  use  Equation  (2)  to  obtain  the 
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equations  for 


3z  3r 


H2  °s2 

be  integrated  from  the  body  (£2  =  0)  to  the  shock  in  order 
to  obtain  z,  x,  and  $  as  functions  of  ?2,  and  £3. 


and 


3(j> 

n - 


These  equations  can  now 


With  the  coordinates  constructed  in  this  manner,  £2  is  always 
orthogonal  to  £^  and  £3  (i.e.,  g23  =  g21  =  0),  but  g^  is  not 
necessarily  orthogonal  to  £3«  By  setting  g22  =  the 

shock  distance,  £2  will  be  orthogonal  to  the  shock  surface, 
and  the  shock  will  correspond  to  €2  *  1.  (Note  that  doing 
this  makes  the  coordinate  system  a  function  of  the  solution 
since  the  shock  distance  is  unknown.) 


Figure  3  illustrates  the  resulting  coordinates  away  from  the 
body. 


Once  z^,^,^),  and  are  known 

everywhere,  the  metrics  g^  can  be  obtained  by  differentia 
tion,  i.e.. 


3z  3z  .  3r  3r  .  2  30  30 


i/j 


1/2,3. 


(5) 


Knowing  the  metrics,  the  following  important  quantities  used 
in  the  PNS  equations  can  be  computed.  These  are  the 
determinant  of  the  metric,  g,  which  for  our  system  is 


9  -  9i]922g  33 


g22g13 


and  the  Christoffel  symbols  of  the  second  kind 


where  the  g*^  are  defined  by 


(6) 
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The  steady  parabolized  Navier-Stokes  equations  in  non- 
dimensional  form  (References  1,2)  can  now  be  written.  In 
order  to  derive  the  PNS  equations,  we  assume  for  the  viscous 

terms  -si—  <<-sr—  and  Jl—  and  that  gradients  of  the  coordinate 
3*1  3*2  °*3 

metrics  are  of  0(1)  .  No  approximations  are  made  to  the  in- 
viscid  terms.  As  an  example,  tne  t ^ -momentum  equation  becomes 
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where  the  velocity  vector  is  written  as 

V  -  U  £x  +  V  a2  +  W  £3 


(10) 
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The  boundary  conditions  at  the  body  are  the  no-slip  conditions# 
a  specified  enthalpy  and  the  pressure  solved  by  using  the 
v-momentum  equation.  The  Rankine-Hugoniot  jump  conditions  are 
used  at  the  shock. 

The  equations  are  solved  by  implicit  differencing  in  the  ^2~^3 
plane.  The  resulting  algebraic  equations  are  solved  by  the 
Newton-Raphson/Gauss  Seidel  iteration  presented  in  References 
1  and  2.  A  plane  of  initial  condition  data  downstream  of  the 
sonic  plane  is  required  in  order  to  obtain  a  solution.  This 
is  obtained  in  the  present  case  by  utilizing  an  asixsymmetric 
shock  layer  solution  for  the  spherical  nose  region  and  trans¬ 
forming  this  solution  into  the  coordinate  system  of  the  PNS 
equations  at  the  sphere-cone  juncture. 

3.  CONCLUSIONS  AND  RESULTS 

The  heat  transfer,  force  and  moment  coefficients  and  surface 
pressure  results  of  four  computed  cases  (Table  1)  have  been 
compared  against  wind  tunnel  data  obtained  by  Arnold  Engineer¬ 
ing  and  Development  Center  (AEDC) .  Selected  results  for  the 
14/7  biconic  with  cuts  and  slice  at  10  angle  of  attack  (Case 
4)  are  presented  here.  Figure  4  describes  the  body  geometry 
used.  The  body  has  a  plane  of  symmetry  and  a  nose  radius  of 
1/2  inch.  Figure  5  compares  the  measured  and  calculated 
force  and  moment  coefficients  for  a  similar  body  but  with  only 
a  windward  cut.  Excellent  agreement  between  the  predictions 
and  the  measurements  are  obtained. 

The  remainder  of  the  comparisons  will  be  surface  pressure  and 
heat  transfer  results  and  will  be  for  the  exact  geomety  of 
Figure  4.  Figure  6  compares  the  heat  transfer  over  the  biconic 
and  cut  regions.  The  axial  distance  is  measured  in  nose  radii 


I 


from  the  apex  of  the  14°  cone.  On  the  windward  side  and  at 
t  90  degrees  the  agreement  is  very  good.  On  the  leeward  side 

up  to  just  past  the  biconic  junction  the  agreement  is  also 
very  good,  but  then  the  solution  starts  to  diverge  from  the 
data.  It  is  felt  that  this  is  due  to  transition  from  laminar 
}  to  turbulent  flow.  For  the  a =10°  case,  separation  in  the 

circumferential  direction  first  occurs  on  the  14°  fore-cone 
near  the  leeside,  but  presents  no  running  problems  for  the 
code. 

> 

Figure  7  is  a  comparison  of  the  wall  pressure  from  the  HYTAC 
calculations  and  the  AEDC  data  and  also  an  inviscid  calcula¬ 
tion  .  The  windward  side  pressures  agree  quite  well,  and 
i  on  the  leeward  side  the  data  and  HYTAC  results  are  also  in 

agreement.  The  axial  distance  in  this  plot  is  in  inches  from 
the  tip  of  the  blunt  body.  Figures  8  and  9  are  circumferential 
plots  of  the  wall  pressure  and  heat  transfer  comparing  the 
f  HYTAC  results  with  the  data  in  the  cut  region  at  an  axial 

distance  of  Z  =  21.77"  (from  the  blunt  nose) .  In  these  plots 
the  X-axis  is  in  inches  from  the  plane  of  symmetry  of  the  body. 

« 

*  The  new  non-orthogonal  body-normal  shock-normal  coordinate 

system  works  extremely  well.  Only  25  points  in  the  normal 
direction  are  needed  to  obtain  good  resolution.  With  19 
circumferential  points  the  solution  over  an  entire  body  is 
^  obtained  with  90  axial  steps  in  approximately  5  minutes  on  a 

CDC  7600.  The  fast  running  time  and  good  agreement  with  data 
is  felt  to  be  due  in  large  part  to  the  new  coordinate  system. 
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Figure  3.  Body  Normal  -  Shock  Normal  Coordinates 
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Force  and  Moment  Coefficient  Comparisons  for  a  14/7  Biconic  with  Windward  Cut 


Figure  7.  I  (AXIAL  INCHES  FROM  NOSE 


Surface  Pressure  Results  Cut  Region  (Z  *  21.77"), 


Figure  9.  Heat  Transfer  Results  Cut  Region  (Z  ■  21.77"), 


